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Abstract 

Real-world sounds always involve partials (or overtones). In this work we describe a mathematic-computational 
method, based on an adaptation of Helmholtz’s theory of consonance, capable of providing an estimation of the self¬ 
consonance and self-dissonance of anharmonic sounds. This method considers, in addition to the given sound to be 
analyzed, a reference harmonic signal with the same fundamental frequency. Consonance and dissonance levels are 
then evaluated for every harmonic component along the reference harmonic sound, also taking into account the partial 
amplitudes, and respective overall indices are obtained. The methodology is illustrated with respect to nonlinearly 
amplified frequency-shifted sounds for increasing levels of anharmonicity, as well as regarding changes in pairwise 
consonance scale relationships implied by anharmonicity. 


‘Human music is the Concordance of divers elements in one 
compound, by which the spiritual nature is joined with the 
body...’ 

Micrologus , Andreas Ornithoparcus (Dowland’s transl.) 

1 Introduction 

Sound and music are inextricably interrelated. While 
the latter can be understood as combinations of sounds 
simultaneous and sequentially, the former consists of a 
seemingly simpler concept. Indeed, at least at first sight, 
sounds seem to be well-defined by their fundamental fre¬ 
quency, involved partials , overall volume , and envelope. 

Two of these properties stand out as being particularly 
complex: partial content and envelope. Interestingly, they 
can be understood as being kind of dual complements, 
because partials have to do with frequency of vibrations, 
while envelope is related to the unfolding of the amplitude 
of these vibrations along time. 

In nature, it is virtually impossible to find a ‘pure’ 
sound, in the sense of involving only its fundamental fre¬ 
quency, without any accompanying partial/overtone. In¬ 
terestingly, also in physics and electronics it is virtually 
impossible to produce such pure sounds. The point is that 
every real-world system involves some unavoidable degree 
of nonlinearity. 

Indeed, a linear system is characterized by its inabil¬ 
ity to produce new frequencies other than those already 


present in the input signal. Contrariwise, non-linearities 
intermix the frequencies in the original signal, producing 
an even more complex partial structure. Similarly, non- 
linearities in sound producing systems also contribute to 
making the partials much more complex not only by com¬ 
bining the original partials, but also by sliding frequencies 
up and down (e.g. with respect to the harmonic series). 

All these effects ultimately converge to making real- 
world sounds very complex regarding their partial content 
and interrelationship. Yet, it is these very same partial 
properties that greatly contribute to defining our resulting 
perception of sounds. 

Human perception is often relative and subjective, and 
may vary along time and space. Yet, there is a common 
set of shared features, that can be understood as defining 
the subsidies for human shared appreciation of sounds. 
In the case of two or more simultaneous sounds, their 
combination can be characterized in terms of the concepts 
of consonance and dissonance (e.g. mm lania ©00). 

As a matter of fact, almost all musical cultures and 
traditions, at least from ancient Indian music to modern 
western approaches have taken these two properties of si¬ 
multaneous sounds into account. In the Bharata classical 
Indian music (e.g. 0 ), for instance, alternating conso¬ 
nance/dissonance is an important subsidy for music com¬ 
position. For instance, it is employed along a musical 
piece in order to enhance specific moods and effects. In 
tonal music, consonance provides the basis for harmony 
and counterpoint (e.g. muni), the two theoretical frame- 
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works dedicated to the study of simultaneous and sequen¬ 
tial, respectively, combinations of notes. 

It is interesting to observe that consonance and dis¬ 
sonance are by no means intrinsically to be sought or 
avoided. In fact, the appreciation of consonance and dis¬ 
sonance usually takes place with respect to some musical 
intention, such as generating surprise and diversity. For 
instance, the leading-note in tonal music (e.g. Hanoi) is 
characterized by strong dissonance while being essential 
for priming modulations and generating ‘expectation’ in 
respective musical pieces. 

Though being human-related perceptual concepts, and 
therefore potentially relative and subjective, consonance 
and dissonance can be approached in a more objective way 
by considering Helmholtz’s respective theory (e.g. [HI 13 

EUSUSI). 

According to this theory, consonance/dissonance stems 
from the combinations between the partials of the two 
involved sounds, especially in which regards beats. Re¬ 
cently [7j, a simple mathematic-computational model 
based on Helmholtz’s theory was reported capable of pro¬ 
viding a surprising agreement with common opinion about 
some of the main temperaments and basic harmony. This 
approach considers not only matching levels between the 
involved partials, but also their respective intensities. 

While consonance approaches can be used to infer pair¬ 
wise combinations of notes, what could them contribute 
to characterizing, through mathematical quantifications, 
the tonal qualities of a single note or sound? 

Harmonic components have been successfully applied to 
the characterization of anharmonicity (e.g. Banana). 
Here, we consider the possibility of extending the concept 
of consonance/dissonance to characterize single notes. 
This approach is aimed at providing information about 
how a single sound can be understood in terms of conso¬ 
nance and dissonance, important concepts in music and 
acoustics. 

A key involved issue regards, given a specific sound, 
how to find a standard reference agains which the partials 
of the given sound can be compared with. A possibility is 
to take some specific model of partials, such as the har¬ 
monic series with the same fundamental as a reference. 

In this work, we specifically address the latter possi¬ 
bility. The hypothesis here is that the harmonic series 
provide a kind of prototype distribution of partials that is 
characteristic of several real-world producing and repro¬ 
ducing sound systems. 

More specifically, it could be posited that humans take 
the harmonic series having the same fundamental fre¬ 
quency as the given sound as a reference for inferring 
its tonal properties in terms of respective ‘consonances’ 
and ‘dissonances’. The harmonic series reference tone is 
henceforth called harmonic reference. The adopted hy¬ 


pothesis is, therefore, directly related to the anharmonic¬ 
ity concept from classical mechanics and molecular spec¬ 
troscopy (e.g mm) as well as Inharmonicity from music 
theory (e.g. Hi HU EJ). 

This article starts by revising the concept of har- 
monicity, and follows by discussing anharmonicity, non¬ 
linear transformations, and frequency-shifted anharmonic 
signals. Then, we present the simple suggested self¬ 
consonance method based on Helmholtz’s respective the¬ 
ory. The potential of this approach is illustrated with re¬ 
spect to the quantification of self-consonance/dissonance 
of progressively anharmonic frequency-shifted sounds 
as well as regarding effects on changing the conso¬ 
nance/dissonance patterns in respective equal tempera¬ 
ment scales. 

2 Harmonic Sounds 

The harmonic sequence defined by the wavelengths lead¬ 
ing to complete periods in a string is given as 

£. — I I I I — (i) 

* 1 ’ 2’ 3’ ’ n’ TV U 

The respective frequencies correspond to 

1,2,3,..., h n ,..., TV (2) 

Where 1 is the fundamental frequency and the remain¬ 
der frequencies are commnly known as overtones or par¬ 
tials. 

Sounds produced by ideal strings (linear dynamics) 
have respective partials obeying the harmonic sequence, 
if /i is the fundamental, the complete sound will have 
frequencies given as 

hifi, h 2 fi, h 3 f i, • • •, h n f x , • • •, h N fi (3) 

Henceforth, any sound (signal) with fundamental /i, 
so that every constituent frequencies h n (i.e. fundamen¬ 
tal plus overtones) belong to the sequence in Equation |3j 
will be called an harmonic signal. We can classify these 
sounds (signals) into two main groups: (i) containing the 
whole set of harmonics, which we will call complete har¬ 
monic sounds (signals); and (ii) containing only a subset 
of the harmonic sequence, which will be called incomplete 
harmonic sounds (signals). 

So, every complete or incomplete harmonic sound with 
fundamental /i will have each of its frequency components 
/ m (assumed to be ordered increasingly) obeying 

fm = K /i (4) 

for some n E {1,2,...} 

However, only complete harmonic sounds will satisfy 
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fn = K /i 


(5) 

for every n G {1, 2,...}. 

For instance, a pipe open at both ends will produce 
(in ideal conditions) complete harmonic sounds, while a 
pipe closed ar one end will generate incomplete harmonic 
sounds respective only to n odd. Sounds emanating from 
ideal strings are complete harmonic sounds. 

3 Anharmonicity and Inhar- 
monicity 

Sounds/signals with overtone structure not fully adhering 
to the harmonic sequence, which is the ubiquitous case 
in the real-world, are henceforth called anharmonic in 
an allusion to the respective classical mechanics concept 
(e.g. p3j). 

As hinted in the introduction section of this work, an¬ 
harmonicity is related to non-linear behavior of dynamical 
systems. For instance, let’s consider a pendulum oscillat¬ 
ing freely from a pivot. The respective dynamics can be 
modeled, according to Newton’s second law of movement, 
as 

6{t) = —a sin(0(t)) (6) 

where 6{t) is the angular position of the bob and a is a 
real-valued constant dependent on the involved physical 
parameters (e.g. [H]). Observe that this ordinary differ¬ 
ential equation is non-linear as a consequence of th sin() 
function. 

For small angular displacements , the above equation 
can be simplified as 

6{t) = —aO(t) (7) 

which is linear and has as solution 6{t) = Asin(ujt + 6) 
for some real parameters A, cj, and e (e.g PH). Observe 
that this type of angular movement 6{t) involves only the 
fundamental frequency, without any overtone. 

Let’s now consider the case of large amplitude oscilla¬ 
tions, which happen when we let the pendulum go from 
an angular position far from the vertical. In this case, the 
solution of the respective differential equation [6] is more 
ellaborate (e.g. m\ yielding the following solution: 

6{t) = 2 asin {s sn [ Ks 2 — ui^t; s 2 ] } (8) 

where sn(x;m), is the Jacobi’s elliptic sine function, Oo 
is the initial position (therefore defining the amplitude of 
the oscillation) and s = sin (Oq/2). 

Figure [l] shows the solutions of the non-linear 
pendulum for several oscillation amplitudes deter¬ 


mined by increasing respective initial positions Oq = 
0.757r, 0.857T, 0.957T, 0.997r. 



Figure 1: Angular displacements 9{t) of a nonlinear pendulum un¬ 
dergoing wide amplitude oscillations corresponding to initial posi¬ 
tions Oo = 0.757T, 0.857T, 0.957T, 0.997T,. 

The resulting oscillations progressively depart from the 
sinusoidal function, acquiring a pronounced non-linearity 
as the amplitude of the oscillation is increased. These 
non-linearities imply several new frequency components 
in the resulting signal, which are not harmonic partials, 
hence the name anharmonic signal. Observe that the ‘fun¬ 
damental’ frequency (corresponding to the largest period 
of the signal) also changes, decreasing with the amplitude 
oscillation. 

The frequency shifts are often proportional to the 
square of the respective amplitude (e.g. [13]), i.e. 

A/ ex A 2 (9) 

In music theory, anharmonicity is alternatively called 
inharmonicity (e.g. pa EH HU). It consists in an im¬ 
portant issue, because every real-world instrument will 
present partials that are not perfectly harmonic. Inhar¬ 
monicity studies this departure from harmonic structure 
unavoidably found in real-world sounds. For instance, pi¬ 
ano inharmonicity has received great attention, as every 
note has ‘stretched’ partials (e.g. PUCES]). The inhar- 
monicities in the piano strings contribute to defining the 
characteristic ‘piano sound’, and therefore needs to be 
taken into account by piano makers and tuners. 

4 Nonlinear Transformations 

Frequently, the music and sounds we hear are the result of 
some electronic amplification , which is required in order 


3 









to produce electric signals with enough power to induce 
audible vibrations in loudspeakers or headphones. 

Even if not always realized, electronic amplification in¬ 
trinsically implies some unavoidable level of nonlinear dis¬ 
tortions. These nonlinearities act on the partials of the 
input sounds, producing many new frequencies. So, a 
particularly interesting way for obtaining altered sounds 
that can be used for consonance/dissonance quantifica¬ 
tion involves producing some input sounds with harmonic 
partials and passing them through a non-linear amplifier. 

We consider this approach in this work. For simplicity’s 
sake, we adopt a quadratic nonlinearity of the type: 



f (Hz) 

Figure 3: The square root of the power spectrum P(f) of the 
squared previous complete harmonic sound x(t ), as obtained nu¬ 
merically. 


y(t) = [Ax(t) + b] 2 (10) 

where x(t) is the input signal, y{t) is the respectively 
amplified output, and A E R is the amplification factor. 
The constant b represents a bias often applied to input sig¬ 
nals in some types of amplification (e.g. class A j20l[2T] h 
This constant determines the operation point of the am¬ 
plification and can be used to control the level of non¬ 
linearity implied by the quadratic amplification. Observe 
that the higher the value of 6, the smaller the nonlinearity 
will be. 

Let’s consider as input complete harmonic signals 
x(t) defined by a respective fundamental frequency /i 
and N partials following the respective harmonic series, 
i.e. /i, 2/i, 3/i,..., n/i,..., N f\. The respective power 
amplitude ai of each of these partials is assumed to follow 
a decaying exponential 


a n = e~ an (11) 

where a E R. Figure [2] illustrates the power spectrum 
P(f) of an input signal x(t) considering a = 0.2, /i = 
110 Hz, and N = 20. 



f (Hz) 


be harmonic signals. This follows immediately from the 
fact that the quadratic transformation (i.e. the product 
of the input with itself) generates only harmonic partials 
corresponding to the sums and subtractions of the original 
harmonics, and other higher powers can be associatively 
understood as successive pairwise products inducing sim¬ 
ilar effects. 

Interestingly, any other nonlinear transformation that 
can be well approximated by a power series (more for¬ 
mally entire analytical functions) will have the same prop¬ 
erty of preserving the harmonicity of the input signal. 
Figure [4] illustrates this fact with respect to the nonlinear 
transformation 


x x 

= 1 + Ti + ¥ 
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( 12 ) 


assuming the same input signal x(t) as in the quadratic 
example above. 



f (Hz) 


Figure 4: The square root of the power spectrum P(f) of the expo¬ 
nential of the previous complete harmonic sound x(t), as obtained 
numerically. 


Figure 2: The power spectrum P(f) of a complete harmonic sound 
x(t) considering a = 0.2, fi = 110 Hz, and N = 20. 

A practical way to obtain the partial content of the 
amplified signal is to apply the discrete Fourier transform 
(e.g. [22]). Figure |3] shows the square root of the power 
spectrum of an output signal y(t) obtained from the pre¬ 
vious input signal for A = 2 and 6=1. 

Given an harmonic signal with fundamental frequency 
/ 1 , any of its polynomial transformations will necessarily 


In an informal sense, the harmonic partial distribu¬ 
tion acts as a kind of 4 eigenfrequencies' 1 of these types of 
polynomial-like transformations (except for the amplitude 
changes and creation of new harmonic frequencies). More 
specifically, using the definitions described in Section [2} 
we have that complete/incomplete harmonic signals will 
be mapped into complete/incomplete harmonic signals. 
These types of transformations will henceforth be called 
harmonic maps. 
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It is important to notice that the above characterized 
harmonic maps is, in principle, restricted to two sounds 
with the same fundamental. In the more general case 
of two distinct fundamentals, new partial frequencies can 
be produced that can strongly influence the respective 
consonance/dissonance (e.g. 7 ). 

So, as we are dealing with given and reference sounds 
with the same fundamental /i, and as harmonic nonlinear 
maps do not change the harmonicity of signals, we need 
to resort to alternative ways of producing anharmonic sig¬ 
nals to test our consonance approach. This issue is ad¬ 
dressed in the next section. 


5 Frequency-Shifted Anharmonic 

Signals 

Let’s us now consider nonlinear transformations capable 
of shifting the harmonic partials. A straightforward way 
to do that is to act directly on the harmonic partials h n = 
1, 2,..., TV, so as to shift their position into new partials 
f n . A possible manner to do that is 

fn = hi (13) 

As the partials are no longer harmonic, combinations of 
respective frequencies as implemented by nonlinear trans¬ 
formations can lead to partials that are not harmonic as 
well as being more dispersed in the frequency domain, 
implying in a multiplicity of new components in the re¬ 
spectively obtained sounds. The higher the value of /3, 
the larger the shifts will be. 

The above approach to obtain anharmonic sounds is il¬ 
lustrated by the square root of the power spectrum shown 
in Figure [5| obtained for the quadratic amplification of the 
same previous harmonic signal (shown in Figure [2j with 
/i = 110 Hz) after having its partials shifted by Equa¬ 
tion [13] for (3 = 1.005. 
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Figure 5: The square root of the power spectra of the signal shown 
in Figure [2] with fi = 110 Hz after having its partials shifted by 
Equation [l3] for fi = 1.005. 

Now, partials appear almost everywhere in the fre¬ 
quency domain despite the relatively small value of /3, so 


that the resulting signal has a substantial level of anhar- 
monicity. In the following section we use this type of sig¬ 
nals in order to illustrate the potential of the consonance- 
based approach to single sound characterization. 


6 The Self-Consonance/Dissonance 
Method 

The principle of the mathematic-computational method 
adopted here in order to quantify the self-consonance 
and self-dissonance of a given sound is outlined in Fig¬ 
ure [6] This method consists in an adaptation of a pre¬ 
viously reported approach aimed at estimating the con¬ 
sonance/dissonance between two distinct notes, which is 
strongly based on Helmholtz’s respective theory (e.g. m 

mm)- 



Figure 6: The framework underlying the adopted method for quan¬ 
tification of the ‘consonance’ and ‘dissonance’ of a given sound, 
whose square root spectrum is shown above ( A ). The harmonic 
reference - below, in green (B) - also with the same fundamen¬ 
tal frequency fi has its harmonic components taken, one by one, 
as the center of two intervals with respective frequency widths Ac 
and Ad. Partials of the given signal B that fall within the former 
interval (shown in blue), i.e. Ac, are considered as ‘in tune’, the par¬ 
tials (identified in red) that are within the interval Ad, but which 
are not in the previous interval, are understood as corresponding 
to beats. These amplitude of these partials are also taken into ac¬ 
count in defining the respective consonance/dissonance. See text 
for additional information. 


The sound to have its consonance/dissonance quanti¬ 
fied ( A ) is henceforth identified as s(£), having funda¬ 
mental frequency fi and respective square root power 
spectrum P(f). The first step consists in obtain¬ 
ing the respective harmonic counterpart signal (F>), 
i.e. the sound that has the same fundamental fi but 
whose partials correspond to the full harmonic sequence 
hifi, h 2 fi, hsfi ,..., ..., h N f 1 . A decaying profile 

of harmonic amplitudes is assumed, as this is often the 
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case with reference real-world sounds (e.g. ideal strings), 
accounting for the relatively less intense contribution of 
higher partials on the resulting consonance/dissonance. 

Two frequency intervals, Ac and Ad, are henceforth 
employed in order to identify the partials of s(t) that 
are here respectively understood to contribute to its self¬ 
consonance and self-dissonance. The rationale for this de¬ 
rives from Helmholtz’s consonance theory (e.g. [niEiE]) 
that understands the overall result of consonance and dis¬ 
sonance the matching (tuning) or not (beats) of partials. 

Similarly to the approach reported in [7], here partial 
differences smaller than 2 Hz are be understood as being 
relatively tuned (the partials rarely last longer than that). 
Diferences between 2 Hz and 60 Hz are here understood as 
potential generators of beats. So, we have that Ac = 4 Hz 
and Ad = 30 Hz. 

Now, for each of the harmonic component hi of the 
reference sound, all the partials fj of the sound A that 
are within the intervals Ac and Ad are identified, and 
the following indices are calculated so as to also take the 
harmonic and partial amplitudes Mi and Mj into account 

C, = Y MiMj \ f\hi — fj\ 

J=H/i hi—fj |<Ac 

di= Y MiMj \f\hi — fj\ 

j=>\fihi-fj\<Ad 
\fihi~fj |>Ac 

In this way, each reference harmonic hi will be assigned 
respective consonance (cf) and dissonance (di) indices. 
Alternative ways can be considered for obtaining the in¬ 
dices Ci and such as not considering the frequency 
difference as being relative to each harmonic component 
(e.g. \fihi — fj | /fihi), considering the square root of the 
magnitudes, among many other possibilities. These alter¬ 
natives can taken into account when adapting the conso¬ 
nance method to particular demands, or tuning it with 
specific styles, etc. Analogue comments can be made re¬ 
garding the choice of the involved parameters. 

Overall indices of consonance and dissonance taking 
into account all the N harmonics can be simply defined 
as 

N N 

c = Y c *; D = Y<ii 

i= 1 i= 1 

One of the interesting advantages of adapting 
Helmholtz’s theory to self-consonance is that now each 
sound will be quantified with respect to both its conso¬ 
nance and dissonance. 


7 Analysis of Anharmonic Time- 
Shifted Sounds 

Figure [9] depicts the consonance and dissonance indices C 
and D obtained for the squared frequency shifted signal in 
the previous example considering several increasing values 
of /? (i.e. increasing anharmonicity). 




Figure 7: Overall consonance (C) and dissonance (D) obtained by 
the proposed method with respect to the squared frequency shifted 
signal in the previous example considering several increasing values 
of (3- A coexistence of scales remindful of fractals can be observed 
in both signatures. 

As expected, the increasing levels of anharmonicity im¬ 
plied by larger values of /? result in decreasing consonance. 
An interesting signature is obtained for the dissonance, 
with alternating peaks and valleys. Both curves have 
a rather intricate structure, with coexistence of several 
scales that is remindful of fractal patterns. Also interest¬ 
ing is the fact that every sound (defined by /3) involves 
degrees of both consonance and dissonance. A good agree¬ 
ment between these estimations and the auditive percep¬ 
tion of the author has been observed. 

The consonance and dissonance signatures obtained for 
the same case as before, but now using b = 0 (less linear 
amplification) are shown in Figure [8j 

The consonance now falls much more steadily (observe 
that the y-axes ranges are different in Figure [9] and Fig¬ 
ure |8|, a consequence of reduced amplification linearity. 
The dissonance signature increases steadily to a substan¬ 
tially much higher level than before. 

The effect of the frequency shifting and square amplifi¬ 
cation can also be inferred by visualizing the consonance 
levels obtained for every pair of sound relationships as 
defined by a given temperament. Figure [To^a) shows the 
consonance relationship estimated for every pair of the 
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beta 

Figure 8: Overall consonance (C) and dissonance (D) for the same 
situation as in the previous example, but now using b = 0. 



Figure 9: The pairwise consonance patterns estimate for the consid¬ 
ered square frequency-shifted sounds as defined for the equal tem¬ 
perament scale. The tonic is at 110 Hz. The self-consonances are 
shown to 1/5 of their original values for the sake intelligibility. Con¬ 
sonances are shown in blue and dissonances in red, and the width 
of the edges reflect the intensity of the consonances. 

12 notes of a equal temperament scale starting at 110 Hz. 
The level of consonance is reflected in the edges width. 
The consonance structure obtained for harmonic sounds 
is also shown in Figure [tO^b) for reference. 

The effect of the anharmonicities in changing the oth¬ 
erwise repetitive consonance patterns can be immediately 
observed. Though the expected cycle of fifths, easily dis¬ 
cernible in (b), is much more tenuous for the anharmonic 
sound relationships in (a). 

8 Concluding Remarks 

A methodology for studying the anharmonicity of sounds 
in terms of self-consonance and self-dissonance based on 
Helmholtz’s theory has been described. Given a sound 
to be studied, a fully harmonic signal counterpart with 
the same fundamental frequency is generated and its har¬ 
monic components are used for identifying respective lev¬ 


els of both dissonance and consonance while also taking 
the respective partial amplitudes into account. 

While other approaches to anharmonicity can provide 
interesting insights to this important phenomenon, the 
Helmholtz theory-based method considered here allows 
the anharmonicity to be expressed in terms of estima¬ 
tions of both consonance and dissonance, which are par¬ 
ticularly important in music theory, practice and appreci¬ 
ation. The application potential of the suggested method¬ 
ology has been illustrated with respect to characteriz¬ 
ing frequency-shifted signals and regarding the estima¬ 
tion of the effects of these alterations in changing the 
consonance/dissonance patterns in an equal temperament 
scale. 

More efforts are necessary to more fully evaluate and 
develop the proposed simple methodology, including the 
consideration of other kinds of anharmonic sounds, as well 
as identification of parameter configurations leading to 
reasonable performance given specific applications. 
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